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Guidance notes 


This computer book contains those sections of the chapters in Block C 
which require you to use Mathcad. Each of these chapters contains 
instructions as to when you should first refer to particular material in this 
computer book, so you are advised not to work on the activities here until 
you have reached the appropriate points in the chapters. 


In order to use this computer book, you will need the following Mathcad 
files. 


Chapter Cl 
121C1-01 Differentiation 


Chapter C2 
121C2-01 Integration as a limit of summations (Optional) 


Chapter C3 
121C3-01 Direction fields and solution curves 
121C3-02 Euler’s method 


Instructions for installing these files onto your computer’s hard disk, and 
for opening them, are given in Chapter AO. 


The computer activities for Chapters Cl and C2 require you also to work 
with Mathcad worksheets which you have created yourself. 


Activities based on software vary both in nature and in length. Sometimes 
the instructions for an activity appear only in the computer book; in other 
cases, instructions are given in the computer book and on screen. 
Feedback on an activity is sometimes provided on screen and is sometimes 
given in the computer book. 


For advice on how each computer session fits into the suggested study 
patterns, refer to the Study guides in the chapters. 


Chapter C1, Section 5 
Optimisation with the computer 


The Mathcad activities for Chapter Cl have only one prepared file 
associated with them. This file explains how to differentiate with the 
computer, which is the topic of Subsection 5.1. ‘To solve optimisation 
problems, in Subsection 5.2, you will be given guidance on how to create 
your own Mathcad worksheets. | 


By the end of this section, you should be able to use Mathcad to carry out 
the process of differentiation and to solve optimisation problems, and also 
be a more skilled and independent Mathcad user. 


5.1 Finding derivatives 


In Mathcad, derivatives are found using the = operator. This operator can 
be used both to find a general formula for the derivative of a function f, 
and to find the numerical value of the derivative of f at a particular point. 


For example, the derivative of the function f(x) = =o can be described by 
the formula f’(a) = 2a, which holds for all real values of 2. Mathcad can 
obtain this expression for the derivative by using the < operator and its 
symbolic commands, as described in Activity 5.1. By defining a value for x 
before these calculations, x = 3 say, the value of the derivative at that 
particular point can be found as f’(3) = 2 x 3 = 1.2, as you will see in 


Activity 5.2. 


Activity 5.1 Finding a formula for the derivative 


Open Mathcad file 121C1-01 Differentiation. Page 1 introduces the 
worksheet. Work through page 2, and then carry out Task 1 on page 3. 


Solutions are given on page 29. 


Comment 


© Each part of Task 1 concerns a function that you were asked to 
differentiate in the main text. In the solutions, where the expression 
obtained using Mathcad differs in appearance from that obtained by 
hand, both expressions are given. If you cannot see immediately why 
such a pair of expressions are equivalent, it is a good idea to copy 
down one expression and use algebra to verify that it can be 
rearranged to give the other. 


© Make sure that the variables entered in the two placeholders of the £ 
operator match. For example, if you mistakenly try to evaluate 
symbolically the expression 


d 
— cos(4xr), 


dt 


then you will obtain the answer 0. 


© Brackets can be crucial when entering an expression into the 
right-hand placeholder of the < operator. For example, if the 
expression x° — 6x7 — 152 + 54 is entered without being enclosed in 
brackets, Mathcad will differentiate the x° term but not the others. 


Remember to create your own 
working copy of the file. 


An example of such a 
quotient is 


sin(t) — t? 
et 


These symbolic keywords 
were introduced in 

Chapter AO, file 121A0-05. 
See also A Guide to Mathcad 
for further details. 


Remember that Mathcad 
notes are optional. 


You should still be working 
with Mathcad file 121C1-01. 


Essentially, this involves a 
direct application of the 
definition of derivative, as 
given by equation (1.3) in 
Chapter Cl. See also the 
second Mathcad note below. 


CHAPTER Ci 


You may find that Mathcad will not allow you to enter into the 
placeholder a quotient whose numerator is a sum of terms (it depends 
on how you enter the quotient). This difficulty can be avoided by 
enclosing the quotient in brackets. 


To avoid such problems, it is good practice to include outer brackets 
around expressions that you enter into the right-hand placeholder of 
the < operator. ‘This also helps to make the expression to be 
differentiated clear on the screen. 

© Sometimes the expression for a derivative obtained by Mathcad can be 
‘improved’ by simplifying it. In place of symbolic evaluation (—), 
either of the symbolic keywords ‘simplify’ and ‘factor’ can be applied 
to obtain a derivative. ‘The outcome from these may or may not be the 
same as that from using — , but will sometimes be in a more 
convenient form. However, what Mathcad regards as ‘simpler’ may not 
necessarily seem so to a human observer! 


Mathcad notes 


The expressions e* and exp(x) are equivalent in Mathcad, but the latter 
form is always used in the output of symbolic calculations, irrespective of 
the form used for input. Similarly, the power + appears in output rather 


2 
than the square root sign. 


Activity 5.1 showed how to use use the eo operator and symbolic 
evaluation to obtain an algebraic expression for the derivative. ‘This 
replicates what you might do by hand, but Mathcad can also\be used to 
differentiate functions that would be rather complicated to do by hand. 
We turn next to how Mathcad can be used to find the numerical value for 
the derivative at a particular point. 


Activity 5.2 Evaluating the derivative at a point 
Turn to page 4 of the worksheet, and carry out ‘Task 2. 


Solutions are given further down page 4 of the worksheet. 


Comment 


© Note the comments made towards the bottom of page 4 of the 
worksheet. Once a value has been defined for the differentiation 
variable (x, say), then the operator can be evaluated either 
symbolically (4(...) —) or numerically (<(...) =) to find the 
numerical value of the derivative at that particular value of z. 


These two evaluation methods usually give the same numerical value, 
but Mathcad calculates the two results in different ways. For symbolic 
evaluation, Mathcad finds a general formula for the derivative and 
then evaluates this formula for the particular value of x, whereas for 
numerical evaluation, Mathcad uses a numerical algorithm to find an 
approximate decimal value. 


© If the expression = f(x) is entered, where a definition for the function 
f(x) is provided earlier in the worksheet, then Mathcad can find either 
a symbolic derivative for f(x) or the numerical value of this derivative 
at any specified point, just as before. 
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Mathcad notes 


© When you enter — to evaluate an expression symbolically, or = to 
evaluate it numerically, it doesn’t matter where on the expression the 
blue editing lines are. All that matters is that the expression is 
complete, with every placeholder filled in. These two evaluation 
methods also behave in a similar way if a change is made to the 
worksheet, above or to the left of a calculation. In automatic 
calculation mode (the default), the result of a calculation involving 
either — or = is updated automatically, while in manual mode, you 
can press the [F9] key to update the result. 


© When evaluating a derivative numerically (<(...) =), you must define 
earlier in the worksheet the point at which the derivative is to be 
found; for example, x := 3. Mathcad then uses a numerical algorithm 
to obtain an approximation to the exact value of the derivative at that 
point, which is usually accurate to 7 or 8 significant figures. Very 
occasionally the method fails, in which case the derivative is 
highlighted in red. Clicking on this expression reveals the error 
message ‘Can’t converge to a solution.’. 


Now close file 121C1-01. 


Sometimes you may want to use Mathcad to find a second-order 
derivative. One way to do this is by using the = operator twice. You can 
enter the = operator in your worksheet, then enter the = operator again 
into the right-hand placeholder, and then fill in all the placeholders 


appropriately. For example, Mathcad gives the result 


+ (Fee — 6x? — 15a” + 54) — 62 — 12. 
Another way to find a second-order derivative in Mathcad is to use the a 
operator, for which there is a button on the ‘Calculus’ toolbar. You should 
enter ‘2’ in the bottom placeholder (which will cause a 2 to appear in the 
top placeholder as well), and fill in the other placeholders just as for the 
operator. You can evaluate the resulting expression either symbolically or 
numerically, in the same way as for expressions involving the 7 operator. 
For example, Mathcad gives the result 

d* 3 2 

aaa (t — 62° —154%+54) — 62-12. 


To change from automatic to 
manual calculation mode, or 
vice versa, click on 
Automatic Calculation in 
the Math menu. 


Without this prior definition, 
x appears in red in the aa 
operator, as an undefined 


variable. 


The keyboard alternative is 
[Ctrl]? (given by the three 
keys [Ctrl], [Shift] and /). 


See Chapter Cl, 
Subsection 2.3. 


See Computer Book A, 
Chapter A3, Subsection 5.5. 


To open a toolbar, either click 
on the appropriate button on 
the ‘Math’ toolbar, or use the 
View menu, Toolbars. 
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5.2 Optimisation 


Optimisation involves finding the greatest or least value taken by a 
function on an interval. In the main text you saw how to apply the 
Optimisation Procedure by hand, but all the calculations required can also 
be performed using Mathcad. In some cases, Mathcad can be applied to 
check calculations already done by hand, while in other cases, calculations 
can be carried out that are too complicated to do by hand. 


In this subsection you will solve two optimisation problems, the first a 
minimisation problem and the second a maximisation problem. 


The orienteer’s problem 


The orienteer’s problem was described in Computer Book A. The orienteer 
starts at a point in a forest and needs to reach a final point on a path. The 
problem is to find where the orienteer should aim to join the path, where 
the joining point is x km from a fixed point O on the path, in order to 
minimise the time taken overall. For the particular data given, this 
involves finding the value of x that gives the least value of the journey time 
(in hours) 


f(x) =0.125V1+4 x? + 0.0625(2 — x) (a in (0, 2)). 
The graph of y = f(x) is shown in Figure 5.1. 


Figure 5.1 Graph of the function for the orienteer’s problem 


In Computer Book A, you found an approximate solution to this problem 
using Mathcad, by zooming in on the graph and applying the graph trace 
tool. You are now in a position to solve the problem more accurately, 
using differentiation. 


In Activity 5.3 you are asked to use the Optimisation Procedure to solve 
the orienteer’s problem. In place of a prepared Mathcad file for this 
activity, you are guided through creating your own Mathcad worksheet. 
This includes a requirement to enter some explanatory text, to make the 
worksheet more comprehensible. 


The instructions for this activity may look lengthy, but several of them 
describe general Mathcad techniques that can save you time and effort, 
now and in the future. The instructions are given mostly in a mouse/click 
way, but details of keyboard alternatives are provided throughout. If you 
wish to use the toolbar buttons, then the activity makes use of buttons on 
the ‘Calculator’, ‘Calculus’, ‘Boolean’ and ‘Symbolic’ toolbars. 


SECTION 5 Optimisation with the computer — 


Activity 5.3 Minimising the orienteer’s journey time 


The following instructions guide you through creating a Mathcad 
worksheet to carry out the Optimisation Procedure, in order to solve the 
orienteer’s problem. Part (b) describes how to enter a title in the 
worksheet. You should similarly enter any other text that you think is 
appropriate. For example, you could include a line of text before each step 
of the Optimisation Procedure, to explain what you are doing. 


(a) 


Begin by creating a new worksheet, as follows. Select the File menu 
and choose New.... In the list of templates that appears, Normal 
should be selected by default. If not, click on it. Then click on the 
OK button to create a new (Normal) worksheet. (Alternatively, type 
[Ctrl]n, or click on the ‘New’ button on the standard toolbar-) 


Enter a title at the top of your worksheet. To do this, click to position 
the red cross cursor in an appropriate place, and choose Text Region 
from the Insert menu. (Alternatively, type a double-quote ", given by 
[Shift]2.) Then type a suitable title, for example, Optimisation - 
the orienteer’s problem, in the text box. To finish, click anywhere 
outside the text box or press [Ctrl] [Shift] [Enter] . If you need to 
edit the text later, simply click on it. 


Enter a definition of the function to be optimised, which is 


f(x) = 0.125V1 + x? + 0.0625(2 — x). 


You can use the keyboard and buttons on the ‘Calculator’ toolbar to 
enter this, or just type the key sequence 


f (x) :0.125*\1+x/A2[Space] [Space] [Space] [Space] +0.0625* (2-x) 
The \ (backslash) keystroke creates the square root sign. 


Now you are ready to solve the orienteer’s problem, by carrying out 
the three steps of the Optimisation Procedure, as follows. 


Step 1: Find the stationary points of f 


Using the = button on the ‘Calculus’ toolbar, enter the expression 
< f(x). Then click on the ‘Equal to’ = button on the ‘Boolean’ 
toolbar, and enter 0 (zero). (Alternatively, just type 

?x [Tab] f(x) [Ctrl1]=0. The ? (question mark, given by [Shift] / ) 
enters the £ operator, [Tab] moves between the placeholders, and 
[Ctr1]= gives the special equals sign. ) 


The stationary points are the points where the derivative of the 
function f is zero, that is, the solutions of the equation just typed in. 
Hence (with the vertical blue editing line still at the right-hand end of 
this equation), click on the ‘solve’ button on the ‘Symbolic’ toolbar, 
then type x into the empty placeholder that appears. (Alternatively, 
just type [Ctr1] [Shift] .solve,x.) Finally, click elsewhere on the 
page, or press [Enter] . 


You should see the solution 0.577 350... appear to the right of ‘—’. 
This is the single value of x at which the function f has a stationary 
point. 


A Guide to Mathcad contains 
detailed information on 
creating and editing your own 
worksheets. See also 

Activity 2.3 in Chapter AO. 


If you have just started 
Mathcad running, then there 
is no need to do this, as it 
automatically starts with a 
new (Normal) worksheet. 


Any text that you enter can 
later be moved, or deleted, as 
you wish. 


If you have difficulty in 
following these instructions, 
then you may like to look 
ahead to the first item in the 
Comment below, which shows 
what should eventually 
appear on the screen. 


Note that this button gives 
the special equals sign, which 
equates the expressions on 
either side of it. (It has 
thicker lines than the equals 
sign obtained from the 
‘Evaluate Numerically’ = 
button on the ‘Calculator’ 
toolbar.) 


A Comment item below 
describes in detail how to use 
the copy and paste facilities. 


Remember that the 
orienteer’s problem is a 
minimisation problem. 


Your worksheet should 
contain the calculations 
shown, but may have different 
text. 


40° 


CHAPTER C1 


Step 2: Evaluate f at each of the relevant points 


Use Mathcad to evaluate the original expression f(x) at the two 
endpoints of the interval, x = 0 and x = 2, and at the stationary point 
between them. For example, type f(0)= to evaluate f(0). It is 
sufficient to input the value for the stationary point to three decimal 
places, x = 0.577, or you can use copy and paste to input all 20 digits 
given in the solution if you wish! 


Step 3: Choose the optimum value 


Identify the least of the three function values that you have found. 
This is the minimum value of f(a) within the given interval. 


You could finish at this point, but it is a good idea to record your 
conclusions in the worksheet. For example, you could enter the text 


The: beast value -of..fG). op [0.21 is.7 Gt -=7). 


replacing the question marks with the values that you found. You can 
add more text if you wish, but the Optimisation Procedure for the 
orienteer’s problem is now completed. 


(e) Finally, save your file, by choosing Save As... from the File menu. 
You will need to give it a suitable name that indicates the contents, for 
example, my121C1-orienteer. This will help you to organise and 
identify your files. 


A solution to part (d) is given on page 29. 


Comment 
© The Mathcad worksheet should now look something like this: 


Optimisation - the orienteer’s problem 
Find the distance x that minimises the time f(x) taken. 
2 , 
f(x) = O0.125,1+x + 00625(2 - x) 
Step 1 Find the stationary points off 


o see) = 0 solve,x — 57735026912962576451 
dx 


Step 2? Evaluate f at the endpoints and at the stationary point 


f(0) = 0.25 f(2) = 0.28 f(0.577) = 0.233 


Step 3 Choose the optimum value 


The least value of ffx) on [0,2] is 0.2353 (at x = 0.577). 


SECT. ION 5) Optimisation with the computer 


© It would have been possible to solve the equation £ f(x) = 0 as 
described above but without explicitly entering ‘=0’. When an 
expression does not contain an equals sign, the symbolic keyword 
‘solve’ will find values of the selected variable for which the expression 
is equal to zero. (However, it is clearer to read such a line of the 
worksheet with ‘=0’ included.) Note also that Mathcad omits the 
zero before the decimal point in the solution. 


© If you want to evaluate f(0.577 350...) using all 20 decimal places 
given in the solution for the stationary point, then rather than enter 
the number from scratch, you can take advantage of Mathcad’s copy 


and paste facilities. To do this, click anywhere in the solution (the Copy and paste is a handy 
horizontal blue editing line extends to select the entire number, way of avoiding the need to 
regardless of where the vertical editing line is positioned within it) and retype awkward or lengthy 
choose Copy from the Edit menu. Then type f() , and paste this expressions. 


number into the empty placeholder between the brackets by choosing 
Paste from the Edit menu. Lastly, click or type = to evaluate the 
function for this value. (The keyboard alternatives for copy and paste 
are [Ctrl]c and [Ctrl]v, respectively. ) 


Mathcad notes 


When working symbolically, a decimal point in the input expression 
triggers a decimal result, with up to 20 decimal places! However, when 
working numerically (evaluating an expression using = ), the number of 
decimal places displayed depends on the setting for ‘Number of decimal 
places’ (Format menu, Result..., ‘Number Format’ tab). By default, the 
results of numerical calculations are displayed to three decimal places. 


Now close the Mathcad file that you have created and saved. 


A traffic planning model 


Traffic planners wish to set up a mathematical model to describe how the 
volume flow rate of traffic (that is, the number of vehicles which pass a 
fixed point in a given time) varies with the average velocity of vehicles 
along a single lane of a road. The eventual purpose of this model is to 
advise on how traffic flow can be maximised. 


The planners assume that each vehicle moves at a constant velocity 

v ms‘. On the basis of a subsidiary model and many observations, they 
estimate that, on average, each driver maintains a distance v + 0.02v? 
metres between the front of their vehicle and the back of the vehicle 
immediately ahead. The average length of a vehicle is estimated to be 

5 metres. The model therefore represents the situation as shown in 


Figure 9.2. 
velocity v 
_ race 
v + 0.02v? v + 0.02v? 
5 5 5 


Figure 5.2 Representation of traffic flow 
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A velocity of 35 ms~! is 
about 126 kilometres per hour 
or 78 miles per hour. 


Remember that this is a 
maximisation problem. 
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The distance between the fronts of successive vehicles is 5 + v + 0.02v? 
metres. Each vehicle, travelling at velocity v ms~', covers this distance in 
(5+ v+4+0.02v?)/v seconds. It follows that the volume flow rate f(v) of 
traffic (in vehicles per second past a fixed point) is given by 
Vv 
Wh 3 ape 


This formula is to apply for 0 < v < 35. 


The traffic planners seek the maximum value of f(v) for v in the interval 
(0, 35]. If such a maximum can be found, it can be used to provide an 
advised speed of travel on the road. 


In the next activity you are asked to create a Mathcad worksheet once 
more, to find both the greatest value of the volume flow rate function f(v) 
and the value of v for which this occurs. 


Activity 5.4 Maximising the volume flow rate of traffic 


(a) Create a new (Normal) worksheet, as in Activity 5.3(a). 
(Alternatively, you may prefer to work with the file that you created 
for the orienteer’s problem in Activity 5.3. To do this, open the file, 
make a working copy of it using a different filename, for example, 
my121C1-traffic, then edit the text and expressions in the worksheet.) 


(b) Enter a title at the top of your worksheet, for example, Maximising 
the volume flow rate of traffic. 


(c) Solve the volume flow rate problem by finding the maximum value of 
the function 
v 


v) = >—————_ 

F() 5+v+0.02v? 
on the interval [0,35] and the corresponding value of v. Do this by 
carrying out the three steps of the Optimisation Procedure, just as you 
did in Activity 5.3(d). 

(d) State your answer in a form which is appropriate in the context of the 
traffic planning model. 


Solutions to parts (c) and (d) are given on page 29. 


Comment 


Having defined the function f(v) for the volume flow rate of traffic in your 
worksheet, you may wish to display a graph of it. You can do so by 
defining a suitable graph range, for example, v := 0,0.1..35, and then 
plotting f(v) against v. While the graph confirms the maximum value at 
v ~ 16, it also shows that there is little difference in the value of f(v) 
between v = 10 and v = 20. 


Save your file, for example as my121C1-traffic. Then close this file. 


Chapter C2, Section 5 
Integration with the computer 


There is only one prepared Mathcad file for this section, and that comes 
towards the end and is optional. As in the later computer activities for 
Chapter C1, you will, for the most part, be creating your own worksheets 
here and using Mathcad directly. Subsection 5.1 shows how to find 
indefinite integrals in Mathcad, while Subsection 5.2 covers definite 
integrals. 


5.1 Finding indefinite integrals 


In Mathcad, indefinite integrals are found using the [ operator. Like the - 
operator, the { operator can be used with Mathcad’s symbolic commands, 
to find an algebraic expression for an integral of a given function. 


In this subsection you are invited to find indefinite integrals for a variety of 
functions. The first activity provides an introduction to finding indefinite 
integrals using Mathcad. 


Activity 5.1 How to find indefinite integrals 


In this activity you will use Mathcad to find the indefinite integral of 2°. 


The buttons referred to below are on the ‘Calculus’ and ‘Symbolic’ 

toolbars. If you wish to use these and they are not already visible, then 

either click on the appropriate buttons on the ‘Math’ toolbar, or select the 

View menu, Toolbars and choose Calculus, then repeat and choose 

Symbolic. 

(a) Create a new (Normal) worksheet. If necessary, see Chapter C1, 
Activity 5.3(a) on page 9 of 


Es eh 3 ksheet, either by clicki th 
(b) Enter the { operator in your worksheet, either by clicking on the [ oe wer sok, 


button on the ‘Calculus’ toolbar, or by using the keyboard alternative 
PCEeriyi . Be careful not to confuse the 


i 
button with the | button, 
(c) Enter the expression to be integrated, x’, in the first placeholder after be ape SA oe 


the integral sign. (This expression is called the integrand.) Then enter  gefinite integrals (as you will 
the variable of integration, x, in the placeholder after the ‘da’. see later). 


(d) Click on the — button (‘Symbolic Evaluation’) on the ‘Symbolic’ 
toolbar, or use [Ctrl]. , the keyboard alternative. Then click 
elsewhere on the page, or press [Enter] , to obtain the integral. Check 
that the answer provided by Mathcad is what you expect. 


(e) Now go through the same procedure to evaluate the integral / u> du. 


(f) If you wish to save your work, then select the File menu and use Save 
As... to name and save your worksheet. (It is a good idea to insert a 
title in your worksheet. If you need to create space for this, then do so 
by positioning the red cross cursor at the top of the worksheet and 
pressing [Enter] to insert as many blank lines as required.) 
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If necessary, refer to the 
instructions in 

Activity 5.1(b)—(d). 

If you arrive at an incorrect 
answer to part (c), then see 
the Mathcad notes below. 
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Comment 


© Notice that the { operator in Mathcad gives only an integral of the 
integrand. It does not give the indefinite integral because it does not 
add an arbitrary constant. 


© The outcomes from integrating [ «dx and f{ u* du demonstrate that 
the form of the indefinite integral depends on the nature of the 
function being integrated but not on the choice of symbol for the 
variable of integration. 


Mathcad notes 


© In part (e), it is sufficient simply to edit the integral expression 
{ x? dz, replacing each ‘a’ by a ‘u’. 


© You can also find an integral of a function f by first defining f and 
then evaluating symbolically the expression [{ f(a) dz. 


© The behaviour of the { operator differs somewhat from that of the £ 
operator, used in Chapter Cl. The result of evaluating the { operator 
symbolically is not affected by a value being defined for the integration 
variable (a, say) earlier in the worksheet. Furthermore, the { operator 
cannot be evaluated numerically, since this makes no sense in the 
context of finding an indefinite integral. If you try this, then Mathcad 
responds with the word ‘function’ to the right of the equals sign. 


In each of the remaining activities in this subsection, you can either 
continue using the worksheet from the previous activity, or close that file 
and create a new worksheet. You will need to decide whether you wish to 
save your work. 


In the next activity, you are asked to use Mathcad to find the indefinite 
integrals of two functions that you integrated by hand in the main text, 
and of a third function which you could also integrate by hand. 


Activity 5.2 Finding indefinite integrals 


Use Mathcad to find each of the following indefinite integrals. 
1 = 
(a) i € + ) fey >) (b) / & = 5sin(5y) dg (y.>D) 


(c) ; (a + cos(axz))dx (where a is a non-zero constant) 
Solutions are given on page 29. 


Comment 


© After the ‘+c’ has been added, the Mathcad answers to parts (a) 
and.(b) agree with those obtained in the main text, and that for 
part (c) is identical to the answer found by applying Table 1.1. In 
part (a), recall that exp(3) is an alternative way of writing e°*. 


© Note that, while Mathcad gives answers here that agree with those 
obtained by hand, there is no simple way of entering the additional 
constraints on the variables that accompany the integrals, for example, 
the condition x > 0 in part (a). For this reason, it is necessary to bear 
in mind that the symbolic manipulations performed by Mathcad might 
not be valid in all circumstances, and that the output should be 
interpreted with care. 


SECTION 5 Integration with the computer 


Mathcad notes 


It is safest to enter all products explicitly in Mathcad, by clicking on the 
‘Multiplication’ x button on the ‘Calculator’ toolbar or by typing *. If 
you enter 3x in part (a), or 5sin or 5y in part (b), then Mathcad will 
assume that you intended to enter a product, and will insert the 
multiplication for you. However, Mathcad will not help in this way if you 
enter ax rather than a*x in part (c). In this case, Mathcad assumes that 
you have entered a single variable name, ‘ax’, rather than the product of 
the variables a and «x. 


The next activity contains four integrals which you found by hand in 
Section 2, followed by three further integrals (in parts (e)—(g)) that cannot 
be found by hand simply on the basis of what is given in the chapter. 


Activity 5.3 Further indefinite integrals 


Use Mathcad to find each of the following indefinite integrals. 


(a) [(e-3)@- de (b) [Ape 


(c) [ sin? da (a) | ade 


ge? +1 


(e) / petit <--@) / oe / mgm tt 


Solutions are given on page 30. 


Comment 


© The answers obtained from Mathcad in parts (a)—(d) are equivalent to 
the answers found by hand. However, they are not always given in an 
identical form, and in part (c) some algebra is required to show the 
equivalence of the two forms. 


© It is possible to use the = operator symbolically to differentiate each 
of the integrals obtained in this activity. This leads to an expression 
equivalent to the original integrand in each case, as expected. 
However, the output obtained is not always identical in form to that of 
the original integrand, and again some algebra is sometimes required 
to show the equivalence of the two forms. 


© Recall (from the context of obtaining derivatives symbolically) that 
the symbolic keywords ‘simplify’ and ‘factor’ can be used as 
alternatives to —. When applied to integrals, either of these may 
again provide an outcome in a more convenient form. 


Parts (e)-(g) of Activity 5.3 show that Mathcad can extend your 
‘integration reach’ beyond the types of functions that you have so far 
learnt how to integrate by hand. However, Mathcad cannot integrate every 
function. For example, if you ask Mathcad to evaluate symbolically either 
of the indefinite integrals 


/ ct a and / Vsin® x da, 


the response is for Mathcad to repeat the given integral without alteration 
(except for the replacement of e~* by exp(—t*) in the first case, and the 
square root by power . in the second). This is how Mathcad responds 
when it cannot find an integral. 


Note that sin? x should be 
input as sin(x)*; for example, 
type sin(x)A2. 


See the final Comment item 
for Chapter Cl, Activity 5.1 
on page 6 of this computer 
book. 


No expression in terms of 
standard functions is known 
for either of these integrals. 
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The [Tab] key provides a 
good way of moving around 
the placeholders. 
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5.2 Definite integrals, areas and summations 


You will need Mathcad file 121C2-01 for (optional) Activity 5.8, later in 
this subsection. First, however, you are invited to create your own 
Mathcad worksheets as in the previous subsection, but now to find definite 
rather than indefinite integrals. 


In Mathcad definite integrals are found using the i operator. This 
operator can be used either symbolically or numerically. If the operator is 
used symbolically, then Mathcad finds an algebraic expression for an 
integral, evaluates this expression at the upper and lower limits of 
integration, and subtracts the second value from the first to find the 
answer. This is the same as the usual approach to finding definite integrals 
by hand. If the operator is used numerically, then Mathcad does not find 
an algebraic expression, but instead uses a numerical algorithm to find an 
approximate value for the definite integral. 


In each of Activities 5.4—5.7 below, as in Subsection 5.1, you can either 
continue using the worksheet from the activity before, or close that file and 
create a new worksheet. You will need to decide whether you wish to save 
your work. 


Activity 5.4 introduces you to the symbolic use of the [ : operator. 


a 


Activity 5.4 How to evaluate definite integrals 


In this activity you will use Mathcad to evaluate the definite integral 


3 
1 
[ <a. 
oO 


(a) Enter the i operator in your worksheet, either by clicking on the _ 


button on the ‘Calculus’ toolbar, or by using the keyboard alternative 
& (the ampersand sign, given on the keyboard by [Shift]7 ). 


(b) Enter the integrand, the variable of integration, and the upper and 
lower limits of integration in the appropriate placeholders. 


(c) Click on the — button (‘Symbolic Evaluation’) on the ‘Symbolic’ 
toolbar, or use [Ctrl]. , the keyboard alternative. Then click 
elsewhere on the page, or press [Enter]. You should obtain the 
answer 


In(3) — In(2). 


Click anywhere on this answer, and then enter =, either by clicking on 
the = button on the ‘Calculator’ toolbar or by typing =, to evaluate it 
numerically. You should obtain the answer 0.405, which is the value of 
In 3 — In2 to three decimal places. 


= 


Comment 


Evaluating symbolically an expression involving the i operator gives an 
expression which is an exact answer (unless the original expression 
contains a decimal point; see the second Mathcad note below). In the 
example in this activity the expression is In(3) — In(2). You can display the 
decimal value of such an expression by evaluating it numerically. This is 
done by selecting the expression, and then entering =. 


SECTION 5 Integration with the computer 


Mathcad notes 


© When you evaluate numerically an expression in Mathcad, the number 
of decimal places displayed is determined by the value of ‘Number of 
decimal places’. The default value of this is 3, but you can change it 
by choosing Result... from the Format menu and then the ‘Number 
Format’ tab. 


© Ifa Mathcad expression involving the £ operator has a decimal point 
in any constant in the integrand, or in both limits of integration, then 
evaluating the expression symbolically gives a decimal answer with up 
to 20 decimal places. (Such an answer is unaffected by the value of 
‘Number of decimal places’.) For example, evaluating symbolically the 
integral below in Mathcad gives the outcome 


see 
/ —dzx — .40546510810816438198. 
o. 


The next activity provides further practice in using Mathcad to evaluate 
definite integrals symbolically. 


Activity 5.5 Evaluating definite integrals 


Each of parts (a)—(d) below gives a definite integral that you were asked to 
evaluate by hand in the main text. In each case, use Mathcad to evaluate 
the definite integral symbolically to obtain an exact answer, and then 
evaluate this answer numerically, to display it as a decimal value. 


= / eee / La ee 
(c) [Sau (d) [ etsineat 


Solutions are given on page 30. 


In the next activity you will see an example of a definite integral that 
Mathcad is unable to evaluate symbolically. When this happens it is often 
worth attempting to evaluate the integral numerically, and the activity 
shows you how to do this. 


Activity 5.6 An awkward definite integral 


(a) Use Mathcad to try to evaluate symbolically the definite integral 


t 3 
/ ee. ah 
0 


You should find that the definite integral is repeated without 
alteration (except for the replacement of e~* by exp(—t?)). This 
means that Mathcad has been unable to calculate an algebraic 
expression for an integral of the integrand, and so it cannot evaluate 
symbolically the given definite integral. 


(b) Evaluate the definite integral numerically, as follows. Click anywere on 
the expression just created, and then enter =. 


You should find that the answer 0.808 is displayed. 


If necessary, refer to the 
instructions in Activity 5.4. 


Remember that 7 can be 
obtained from the ‘Calculator’ 
or ‘Greek’ toolbar, or by 
typing [Ctrl] [ShiftIp. 


At the end of Subsection 5.1 
it was pointed out that 
Mathcad cannot evaluate the 
indefinite integral 


fat dt. 
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Activity 5.8 indicates a 
possible basis for such an 
algorithm. 


The first two steps here are 
identical to those for symbolic 
integration. ‘These steps 
would also have been required 
for the definite integral in 
Activity 5.6(b), had the 
expression not previously 
been entered. 


For example, the values of all 
the definite integrals in 
Activity 5.5 can be found 
satisfactorily using numerical 
integration. 
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Comment 


© To evaluate any definite integral numerically, you should enter it in the 
same way as for symbolic evaluation, then select it and enter =. 


© The reason why some definite integrals can be evaluated numerically 
but not symbolically in Mathcad is that symbolic evaluation requires 
Mathcad to find an algebraic expression for an integral, whereas 
numerical evaluation involves the use of a numerical algorithm. The 
answer obtained from this algorithm is an approximation, though 
usually an accurate one. 


Mathcad notes 


On rare occasions, the numerical method used by Mathcad for evaluating 
definite integrals fails to produce a value. In such a case, the integral is 
marked with the error message ‘Can’t converge to a solution.’. 


The next activity gives you further practice in using the rs operator 
numerically. 


Activity 5.7 Evaluating definite integrals numerically 


1 
(a) Evaluate numerically the definite integral / t? dt, as follows. 
0 


(i) Either click on the ae button on the ‘Calculus’ toolbar, or type & 
to insert the definite integral operator and its placeholders. 


(ii) Fill in the four placeholders appropriately. 


(iii) The blue editing lines should be on the expression; if not, select 
some part of the integral. Then evaluate the expression numerically, 
by entering =. 


1 

(b) Evaluate numerically the definite integral { V1— 2? dz. 
<4 

Solutions are given on page 30. 


Comment 


These two definite integrals can be evaluated either numerically or 
symbolically. As you can check, symbolic evaluation gives = for (a) and 7 
for (b). 


You have seen that in Mathcad many definite integrals can be evaluated 
either symbolically or numerically. You might wonder which it is more 
appropriate to invoke in any given situation. If you want an exact answer 
(so you can see where constants such as 7 feature in it, for example), or if 
a general result is required, such as a formula for 


1 
/ cos(kx)dx (where k is a non-zero constant), 
0 


then use the symbolic approach. If you simply want a number that is an 
accurate value for the definite integral, then numerical evaluation should 
suffice. If you want both an exact value and a decimal value for the 
answer, then you can first evaluate the definite integral symbolically and 
afterwards enter =. 


SECTION 5 Integration with the computer | 


Integration as a limit of summations 


In Subsection 4.2, you saw that a definite integral could be approximated The remainder of this 
as closely as required by a finite sum. This was demonstrated in particular subsection will not be 
for the case in which the value of the definite integral gives the area assessed. 

beneath the graph of a function, and each finite sum represents the overall 

area of a set of rectangles. Each rectangle is based upon a subinterval, and 

the approximation to the definite integral improves as the number of 

subintervals is increased. 


For the area beneath the graph of the function f(x) between x = a and We assume (as in the main 
x = b, which is given exactly by the definite integral text) that f(a) takes only 
5 non-negative values in the 

/ f(a) dx, interval |a, b]. 
the approximation based on N subintervals is 

N21 N-1 ig 

we hf(a+ih)=h ‘S f(a+tih), where h= W See equation (4.1) in 

i=0 i=0 Chapter C2, Subsection 4.2. 


The remaining (optional) activity in this section asks you to use the 
prepared Mathcad file 121C2-01 to explore the relationship between these 
approximations to the area under the curve and the definite integral itself. 


. Activity 5.8 Integration as a limit of summations (Optional) 


(a) Open Mathcad file 121C2-01 Integration as a limit of 
summations, and read through the worksheet, which consists of a 
single page. The definite integral being approximated here is 


40 
/ 15,/sin (= dx, 


whose value was sought by a process of successive approximation in 
Subsection 4.2. 


The value of this integral is the area on the graph beneath the red 
curve. The value calculated for A is an estimate for this area, using 
approximation by rectangles based on N subintervals. The value of NV 
is initially set to 4. 


Investigate the effect of increasing the number of subintervals, NV. Use 
in turn the following values for N: 


20, 50, 100, 500, L000, 5000, 10000. 


Compare the corresponding values obtained for A with those in the 
right-hand column of the table below. 


Number of Sum of areas 
subintervals of rectangles 
4 402.27 
20 42.1 This is Table 4.1 in 
50 A56.41 Chapter C2, Subsection 4.2. 
100 AN 721. 
500 A5T.62 
1 000 457.64 
5 000 457.66 
10 000 457.66 
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(b) You may like to use the worksheet for other definite integrals, to 
investigate the behaviour of area estimates A as N is increased. For 
example, you could investigate these estimates for the definite integral 


1 
mae 
[ eae, 

0 


by first editing the worksheet so that f(x) =e~* and b=1. 


If you have time, look also at the behaviour of the corresponding area 
estimates for the definite integral 


m/A 
/ tan © dx. 
0 


Comment 


© You should find that the numerical values obtained for A in part (a) 
match to two decimal places the values given in the right-hand column 
of the table. These area estimates appear to tend to the limit 457.66 
(to two decimal places), which is also the numerical value provided by 
Mathcad for the definite integral. 


© The values of the definite integrals suggested in part (b) are 


1 4 m/4 
/ e* G@r=Use0s and / tana ge = UAT, 
0 0 


each to three decimal places. In each case, the estimates A converge to 
this value as N increases. 


In the first example, the area estimates are greater than the value 
given by the definite integral for the area under the curve, while in the 
second example they are smaller. You can see this illustrated on the 
graph by setting small values of N (N < 20, say). In the first example, 
the tops of the rectangles all lie above the curve, while in the second 
example they all lie below. 


Mathcad notes 


© The summation sign is obtained from the = button on the ‘Calculus’ 


toolbar, or by typing [Ctr1]$ (for which you have to press the three 
keys[Ctr1] , [Shift] and 4 together). 


© Note that the definite integral is set up in this file with integrand f(z), 
where the function f(a) is defined earlier in the worksheet. This 
approach is possible for either numerical or symbolic evaluation of a 
definite integral. 


© The rectangles are filled in by drawing zig-zag lines very close together. 
On the screen these lines give the appearance of a solid block of colour, 
but the tiny gaps between the lines may become apparent if printed. 


Now close Mathcad file 121C'2-01. 


Chapter C3, Section 4 
Differential equations with the computer 


In this section you will use the computer to see how the information 
contained in a differential equation can be displayed graphically, and how 
differential equations can be ‘solved’ numerically and graphically, even 
where no formula for the solution can be found. 


There are two prepared Mathcad files associated with this section. ‘The 
first draws direction fields and plots solution curves (graphs of solutions) 
for first-order differential equations, and the second illustrates how a 
numerical method for obtaining the solution to an initial-value problem 
works in practice. 


This numerical method can be applied to any differential equation of the 


form 
d 
= diz), As special cases, f may be a 
d function of x alone or of y 


where f is a known function of two variables. This form includes the types alone. 
of differential equation considered in the main text, but also others. 


Mathcad does not contain any symbolic facilities for solving differential 

equations directly, to find a formula for the solution. However, as you saw 

in Chapter C2, Mathcad can be used to find integrals, and finding integrals 

is the main constituent of the two methods for solving first-order 

differential equations that are introduced in the main text (direct 

integration and separation of variables). Hence Mathcad can be used 

indirectly to help solve some first-order differential equations. Mathcad 

can also be applied to differentiate a solution that you have already found, Differentiation was the 
to check that it is indeed a solution to the given differential equation. subject of Chapter Cl. 


In this section, however, Mathcad will be used numerically and graphically 
rather than symbolically. 


4.1 Direction fields and solution curves 


The differential equation 


dy 


cannot be solved, using the methods of this chapter, to give an equation 
relating x and y. However, the direction field of this differential equation 
provides enough information to indicate the different types of solution 
curve that occur. In Activity 4.1 you will see this direction field drawn by 
Mathcad, and will be able to ask Mathcad to plot the solution curve 
through any point of your choice. 
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Activity 4.1 The differential equation dy/dx=x+y 


Open Mathcad file 121C3-01 Direction fields and solution curves. 
The worksheet opens with the direction field of the differential equation 
dy 


rR . 


drawn for a grid of points with integer coordinates in the (x, y)-plane, for 
values of x between —5 and 5, and values of y between —2 and 4. For each 
such point (x,y) a line segment is plotted through the point, and the slope 
of this line segment is f(x,y) =x+y. 


(a) Set S to 1, so that a solution curve is plotted through the point (0,0). 

Briefly describe the curve obtained, or make a small sketch of it. Now 

The initial condition is change the value of yo, to obtain a solution curve through each of the 
y(Xo) = Yo. following points in turn: 


(6d), 109 Gea, ie 


In each case, note down a brief description of the curve or make a 
small sketch of it. 


(b) Can you group the solution curves seen in part (a) into distinct types 
of behaviour? How many different types of behaviour are there? 


For each of the following points, try to predict from looking at the 
direction field which type of behaviour the solution curve through the 
point will exhibit: 


(—3,-1), (—1,0), (4,2). 


Then, by making a suitable choice of values for 7 and yo, use 
Mathcad to plot the corresponding solution curve and to confirm your 
prediction. 


Solutions are given on page 30. 


Comment 


The line segments of the direction field give a good indication of where the 
solution curves lie. 


Mathcad notes 


© A definition for a function of two variables is created in the same way 
as that for one variable. For example, to define the function 
f(x,y) =x+y in Mathcad, you could type f(x,y) :xty. 


© The calculations used to draw the direction field and solution curve 
are ‘hidden’ off the page of the Mathcad worksheet to the right. The 
area beyond the right-hand margin of a Mathcad page (which is 
marked by a solid vertical line) can be used just like the rest of the 
worksheet. It is divided into further pages, where you can place 
mathematical expressions, text, graphs and pictures. You do not need 
to look at the calculations in this worksheet, but in general, you can 
view what is in the ‘hidden’ area by using the horizontal scroll bar to 
move to the right. 


In the next activity you are asked to use Mathcad to plot the direction 
fields for two other first-order differential equations. In each case, you can 
try to visualise from the direction field how the solution curves behave, 
before plotting some of these curves. 
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- Activity 4.2 Direction fields and solution curves S 


(a) 


Investigate the direction field and solution curves for the first-order 
differential equation 


dy 

ae COs. 2 = + 

dx ‘ : 
as follows. 


(i) First set the variable S to zero, so that no solution curve is 
plotted. Then enter the right-hand side for this differential equation 
into the definition for f(a, y). 


(ii) Try to predict from the direction field where the solution curves 
lie. What types of behaviour will the solution curves exhibit? 


(iii) Set the variable S to 1 to display a solution curve, and then try 
different values for x) and yo to confirm the predicted behaviour of 
solution curves. Use your observations to describe the behaviour of the 
solution curves for the differential equation. 


Now follow the same procedure as in part (a) to investigate the 
direction field and solution curves for the first-order differential 
equation 

dy 

= a ty, 
However, in this case, start by altering the scope of the grid in the 
y-direction, by setting Yl := —5, Y2:=5 and q:= 10. 


Solutions are given on page 31. 


Comment 


° 


The differential equation in part (a) is of the form dy/dx = f(x), and 
so can be solved in principle by direct integration. However, it turns 
out not to be possible to find an algebraic expression for 


[cm 1a. 


The presence of the cosx means that the slopes of the direction field 
repeat at horizontal intervals of 27. (The slopes are also invariant in 
the vertical direction, for any given choice of x, as noted in the 
solution.) The alternate positive and negative slopes indicate that 
solution curves will undulate. The overall increasing trend is not so 
obvious from the direction field, but might be expected because the 
magnitude (steepness) of the slopes appears to be greater where the 
slopes are positive than where they are negative. 


The differential equation in part (b) cannot be solved by the methods 
of this course. Using Mathcad to plot the direction field provides a 
good approach to visualising the types of solution curve for such a 
differential equation. 


Now close Mathcad file 121C8-01. 


You should still be working 
with Mathcad file 121C3-01. 


This was the case for each of 
the differential equations in 
Activity 4.2. 


Recall that a direction field 
provides a slope value at 
every point within a given 
region, and not just at the 
particular grid points 
displayed on a graph. 


You studied this differential 
equation in Activity 4.2(a). 


In the file, the initial values of 
x and y are denoted 
respectively by xp and yo. 
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4.2 Euler’s method 


In Subsection 4.1 you saw how a direction field can be used to visualise the 
information provided by a first-order differential equation. You also saw 
solution curves plotted on top of the direction field. It is straightforward to 
plot such curves for a first-order differential equation whose solutions can 
be expressed in terms of a simple algebraic formula. However, even where 
no such formula can be found, approximate solution curves for a 
differential equation can still be plotted. These graphs are based on a 
sequence of numerical estimates for solution values, which constitute a 
numerical solution to the differential equation. In this subsection you will 
see how such a numerical solution can be obtained. 


The procedure to be described below, for obtaining an approximate 
numerical solution to a first-order differential equation, is known as Euler’s 
method. To see how the solution is built up step by step, it is illuminating 
to consider the corresponding graphical construction. This involves 
forming a connected chain of line segments, each of which has a gradient 
given by the slope of the direction field at the left-hand end of the line 
segment, as shown in Figure 4.1. 


slope = f(2,y) 


Figure 4.1 Graphical construction for Euler’s method, where dy/dx = f(x,y) 


You will see this idea explained in greater detail in the next activity. 


Activity 4.3 Introducing Euler’s method 


Open Mathcad file 121C3-02 Euler’s method, read the introduction 
and turn to page 2 of the worksheet. Here a direction field has been set up 
for the function f(x,y) = e°** — 1, that is, for the differential equation 


dy 
a ee Ps iz 
dx 
Also, an initial condition is specified, as 


y=OQwhens = thatis, (6) = 0. 


Note, to the right of the graph, that the value f(r, y,) = 1.718 is 
displayed, where n = 0. This is (to 3 d.p.) the value of 


f(0,0) =e° -l=e-—l, 


which is the slope of the direction field at the starting point 

(20, yo) = (0,0). This starting point is denoted on the direction field by a 
small blue box, and the direction of the direction field at this point is 
coloured magenta. 


SECTION 4_ Differential equations with the computer 


You will now see how an approximate solution to this initial-value problem 
can be built up graphically, step by step. Ensure that all of page 2 from 
the heading ‘Solution curve’ to the bottom of the graph is visible on your 
screen. 


(a) Change the value under ‘Number of steps’ in turn to N = 1, N = 2, 
N =3, N =4 and N = 5. In each case, observe the effects that the 
change in value causes. 


(b) Observe the effect of changing, in turn, the step size to h = 0.5 and 
the number of steps to NV = 10. 


(c) Observe the effect of changing the initial values. For example, set 
fy) = 1 ene we = 2. 


Comment 


© When ‘Number of steps’ is changed to N = 1, the first segment of the 
approximate solution curve is drawn, from the starting point 
(Zo, Yo) = (0,0) to the point (x1, y,) = (1, 1.718). The coordinates of 
these points appear in the tables to the right of the graph. The slopes 
of the direction field at these points, f(z, yo) = 1.718 and 
f (x1, y1) = 0.717, are also shown. 


The two broken blue line segments (one coinciding with the z-axis) 
that appear on the graph are temporary construction lines. They 
illustrate how the first segment of the solution curve is drawn as the 
hypotenuse of a right-angled triangle. The base (run) of this triangle is 
equal to the step size, which is specified above the graph as h = 1. 
The triangle is constructed, as shown in Figure 4.2, so that the 
gradient of its hypotenuse is equal to the gradient given by the direction 
field at its left-hand vertex. This gradient is f(2%o, yo). On the other 
hand, the gradient (slope) of the hypotenuse is given as usual by the 
rule ‘slope equals rise over run’, where the run is the step size, h. The 
rise, which is the height of the triangle, is therefore given by 


rise = run X slope = hf (Zo, yo). 
The coordinates (x1, y;) of the top vertex of the triangle can then be 


calculated from the coordinates (29, yo), the rise and the run. This 
gives 


(x1, Y1) = (xo + run, Yo + rise) 


= (xo = a h, Yo + hf (xo, Yo)) 
= (0+ 1, 041.718) = (1, 1.718). 


© When the number of steps is changed to N = 2, the second segment of 
the approximate solution curve is added, from the point 
(v1,y1) = (1, 1.718) to the point (x2, y2) = (2, 2.435). The broken blue 
lines now illustrate the ‘construction triangle’ for this second line 
segment, as shown also in Figure 4.3. The base (run) of this triangle is 
again equal to the step size, h = 1, but the gradient of its hypotenuse 
is equal to the gradient given by the direction field at (21, y,). This 
gradient is f(2,,y,) = 0.717, as indicated by the table on the right of 
the screen. The coordinates (22, yo) are therefore given by 


(Yo, yo) = (%, + run, y; + rise) 
= (0, +h, we + Bf leg) 
S(t 77 Oi — 12 2456). 


Comment on part (a) 
for N =1 


All of the values are shown to 
three decimal places. The 
slope of the direction field at 
(a1, y1) is now shown by a 
line segment on the graph, 
coloured magenta. 


(x1, 41) 


hf (xo, yo ) 


Figure 4.2 First segment 


Comment on part (a) 
orn =F 


Figure 4.3 Second segment 
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Comment on part (a) © 
ion NV ee 3, VV = 4nd Nh SS 


Note that the values of 

f (Xn, Yn) given in the table 
are negative for 

n = 2,3 and 4. These 
correspond to the negative 
slopes of the direction field at 
the corresponding points 

(Xn, Yn), as is clear from the 
graph. 


(Tati, Yu+i) 


hf (xn, Yn) 


Figure 4.4 (n+1)th 


segment 
Comment on part (b) © 
Comment on part (c) © 


With two steps now completed, it should be fairly clear how the 
process continues when the number of steps is changed to N = 3, 

N =4 and N =5. The chain of line segments (which form the 
approximate solution curve) and the rows of corresponding values in 
the tables to the right of the graph build up one by one as N is 
incremented. Each time, the latest line segment is joined to the 
previous one, and its slope matches that of the direction field at its 
left-hand end. 


Mathematically, this means that in constructing the line segment from 
(%.. Uo) tO a), we ave 


(inv, Vai) = (n+ ranyy, + rise) 
= (6, 8, Gn FPF (Gn, 9n)) (ee = 0,1, 2,-.-), 


This is illustrated in Figure 4.4. In other words, the sequences z,, and 
y,, are determined by the pair of recurrence relations 


Un+l = In +h, Yn+1 =F, EPs) be: 07 1)2;-..), 


as shown to the right of the graph on screen, together with the 
specified values for x) and yo. These recurrence relations encapsulate 
Euler’s method. 


Each value y,, is an estimate of the ‘true solution’ y at x = ©; that is, 
y, is an estimate of y(z,). Clearly, the sequence of estimates obtained 
depends on the choice of both the step size, h, and the number of 
steps, NV. 

With N steps of size h, Euler’s method provides N + 1 solution 
estimates, spaced at regular horizontal intervals between the chosen 
starting value, 2), and 79 + NA. 


Note that, when starting at (zo, yo) = (0,0), either 5 steps of size 1 

or 10 steps of size 0.5 provide a final solution estimate at 7 = 5. By 
halving the step size and doubling the number JN of steps, their 
product Nh, and hence also 2 + Nh, remains the same. Notice, 
however, that the solution estimates obtained at x = 5 are different: 
the final coordinates are (x5, y5) = (5, 0.986) in the first case and 
(10,410) = (5, 0.558) in the second. Halving the step size and 
doubling the number of steps improves the accuracy of the solution 
estimate. In the next activity you can investigate further how accuracy 
improves with reductions in step size. 


Setting x) = 1 and yo = 2 produces a different chain of line segments, 
which is again an approximate solution curve for the differential 
equation, now with the initial condition y(1) = 2. (The same is true 
for other choices of initial values.) If h and N are unaltered from 
part (b), then the last row of values at the right of the screen 
corresponds to a point which lies outside the displayed graph region. 


Activity 4.3 showed how Euler’s method works. You saw a connected chain 
of line segments built up, one by one, as an approximation to a solution 
curve. However, this was a fairly inaccurate approximation. The 
construction utilises slope values provided by the direction field only at 
points which are a horizontal distance h apart, where the values assigned 
to the step size h in Activity 4.3 were first 1 and then 0.5. As a result, the 
approximate solution curve was based on very limited information. 
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More information can be extracted from the direction field by reducing the 
step size, provided that also the number N of steps is increased to 
maintain coverage of the x-values over which a solution is sought. The next 
activity shows that such use of extra information leads to improvements in 
accuracy, and that an estimate for the solution y(x) at a particular chosen 
value of x can, in principle, be found to whatever accuracy is required. 


_ Activity 44 Using Euler’s method 


Turn to page 3 of the worksheet, and read the first paragraph. Then scroll You should still be working 
down until all of the page from the heading ‘Graphical and numerical with Mathcad file 121C3-02. 
solution’ to the bottom of the graph is visible on your screen. 


The page is set up to solve the same initial-value problem that was — 
considered in Activity 4.3(a) and (b), namely, 


dy 
=e 4 0) =0. 


The computation, using Euler’s method, takes place from the starting 

value x = Zo to the finishing value x = xval, where xval can be specified 

by the user. ‘The step size h can be chosen as before. However, in contrast 

to the situation in Activity 4.3, it is not possible here to vary 

independently the number of steps, NV. Given a step size h, the number of 

steps to be used is calculated automatically from the condition that the 

final step is to reach x = xval. As a result, the value calculated for yy is In other words, if h divides 


always an estimate for the value y(xval) of the true solution at x= aval. exactly into xval — xo, then 


(a) Change the value of xval to 9 (from its starting value of 5). What N = (xval — xo)/h. 


changes on the screen, and what stays the same? 


(b) We now seek an estimate for the value of y(9) (the true solution value 
at x = 9). This is obtained by progressively reducing the step size. 


First note the value of yy that appears on screen (this is the estimate 
for y(9) obtained with step size h = 1, and hence with 9 steps). 


Now change the step size in turn to 0.5, 0.2, 0.1, 0.01, 0.001 
and 0.0001. In each case, note the corresponding value for yy. Use 
these values to estimate the value of y(9) to one decimal place. 


Solutions are given on page 31. 


Comment 


© As the step size A is decreased, the number of steps increases. With 
h = 0.0001, there are 90000 steps, for which the calculation may take 
an appreciable time on your computer. 


© With h = 0.5, the approximate solution curve still looks like a 
connected chain of line segments. However, with h = 0.2 the graph 
appears significantly more like a smooth curve, and this remains the 
case for smaller values of h. (The graph is still in fact made up of 
short line segments, but so is every graph drawn by Mathcad with the 
trace type set to ‘lines’ !) 


© Just as the numerical estimates for y(xval) appear to converge as the 
step size is decreased, so too do the approximate solution curves seem 
to converge towards a ‘limiting curve’ on the direction field. The 
graphs obtained provide increasingly close representations of the 
actual solution curve. 


Zi 


You should still be working 
with Mathcad file 121C3-02, 
on page 3 of the worksheet. 
Before doing anything else, 
reset the step size toh =1 
(see the first item in the 
Comment below). 


See the Comment for 
Activity 4.4 of Chapter B1, in 
Computer Book B. 
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Euler’s method, as illustrated in Activities 4.3 and 4.4, was used to plot 
the solution curves in file 121C3-01, used for Activities 4.1 and 4.2. There, 
the number of steps was set to N = 1000, and calculations were made to 
draw the line segments to the left, as well as to the right, of the initial 
point. In that way, the line segments gave the appearance of a smooth 
solution curve over the whole horizontal graph range. 


If you are interested and have the time, then you might like to try the 
following optional activity, which involves a function f(x,y) that depends 
on the dependent variable y as well as on the independent variable x. For 
this particular example, it is also possible to check the outcome against a 
formula for the solution. 


Activity 4.5 Another initial-value problem (Optional) 


(a) Use Euler’s method to estimate to three decimal places the value 
of y(6), where y(x) satisfies the initial-value problem 
dy 


ees ret 
x Mes tel al y(1) 


(For the step size h, use in turn the values 1, 0.5, 0.1, 0.01, 0.001 
and 0.0001.) 


(b) Check that the function y = x — 4+ 4e*~* is the solution to the 
initial-value problem in part (a). Hence find y(6) exactly. Does this 
agree with the value that you obtained using Euler’s method in 
part (a)? 


Solutions are given on page 32. 


Comment 


© If you do not start by setting h = 1, then any other change will cause 
recalculation for the most recently-used value of the step size, 
h = 0.0001, which may be time-consuming. Alternatively, recall that 
any Mathcad calculation can be interrupted by pressing [Esc] and 
then clicking ‘OK’ in the resulting option box. You may prefer to 
change here to ‘manual calculation mode’. 


© There is no need to alter any part of the worksheet before page 3, nor 
to change the parameters which define the grid for the direction field. 
On page 3 you need to alter the direction field function definition to 
f(x,y) =2—y—3. Also, the values of zp and yo should both be set 
to 1, and the value of xval to 6. 


© Using the suggested values of h in turn, the numerical estimates yy 
appear to converge and the approximate solution curves do likewise. 
The ‘chain of line segments’ is visible for h = 1 and for h = 0.5, but no 
departures from smoothness are apparent on the graph for smaller step 
SIZeS. 


Now close Mathcad file 121C3-02. 


Solutions to Activities 


Chapter C1 


Solution 5.1 


Where more than one expression is given below for a 
solution, the first is similar to the Mathcad output 
and the second is a form that you are more likely to 
obtain by hand. (You found each of the derivatives 
by hand in the main text of Chapter C1, as indicated 
by the references below.) 
(a) 327-122 —15 
See Activity 2.2(a). 
(b) 4ar? 
See Activity 3.4(a). 
es (cos(t) — 2t) (sin(t) — t?) 
exp(t) exp(t) 
cost — sint +t? — 2¢ 
a 
See Exercise 4.2(b). 
(d) 2cos(x*)z = 22 cos(zx?) 
See pages 48-49. 
(e) —4sin(4x) 
See Activity 4.7(a). 


(t? + 3) 
t 


(f) 2tIn(t) + 
See Activity 4.2(c). 
1 In(u) 
(8) u(u2+3) ~ (u2 +3)? 
See Activity 4.4(b) 
(exp(t)+1) e'+1 


= 1 +3 — Pu? Inu 
u(u? + 3)? 


(h) (exp(t) + t) ve 
See Exercise 4.3(d). 
Solution 5.3 


(d) Step 1: The one stationary point of the function 
f(x) =0.125V/ 1+ x2 + 0.0625(2 — zx), 
with domain [0,2], is at « = 0.577 (to 3 d.p.). 


Step 2: The values of f at the interval endpoints 
are 


f(0) = 0.256. aad ((2) = 0.280, 
while the value of f at the stationary point is 
JOaTT) =O Gate 3 dp... 


Step 3: Hence the minimum value of f(«) for x 
in the interval [0,2] is 0.233 at x = 0.577 (both 
to 3 d.p.). 


Hence the solution to the orienteer’s problem is 
to aim to join the path at a distance 0.577 km 
from the fixed point O on the path. 


(This agrees with the answer 0.58 km found in 
Activity 5.7(a) of Chapter A3 in Computer 
Book A.) 


Solution 5.4 
(c) Step 1: The stationary points of the function 
v 


WO) = ST 000 
ate at v= +12 (ta 3 4 p.); 


Step 2: Only the positive stationary point, 
15.811, lies inside the interval [0,35]. The values 
of f at the interval endpoints are 


f(0)=0 and (35) = 0.543, 


while the value of f at the stationary point in 
the interval is 


f (15.811) = 0.613 (all to 3 d.p.). 


Step 3: Hence the maximum value of f(v) for v 
in the interval [0,35] is 0.613 at v = 15.811 
(both to 3 d.p.). 


(d) According to the model, a maximum traffic flow 
rate of about 0.61 vehicles per second can be 
achieved, by keeping the speed of traffic at 
about 16 ms~! (that is, about 57 km per hour 
or 35 mph). 


Chapter C2 


Solution 5.2 


In each case, c is an arbitrary constant which has 
been added to the expression given by Mathcad. The 
answers to parts (a) and (b) agree with those 
obtained earlier by hand. 


(a) / € + *) dx = In(x) + 4 exp(3z) +¢ 
See Activity 1.2(a). 

(b) / € + 5sin(6y) dy = “5 — cos(5y) + ¢ 
See Exercise 1.1(b). 


sin(ax) 


(c) Jo + cos(ax)) dx = ax + 
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Solution 5.3 


In each case, c is an arbitrary constant which has 


been added to the expression given by Mathcad. The 


Mathcad answers in parts (a), (b) and (d) resemble 
closely those obtained by hand in the main text. 


(a) [(e-3)(@-1)de =} Ly? _ In* + 344+¢ 
See Example 2.1(a). 
203 
Sx 


See Example 2.1(c). 


dg = 4 73/2 _ 6z/? + ¢ 


(b) 


(c) [sw a dx = —}cos(z) sin(x) + 54 +c 


In Activity 2.3(a), the answer $2 — 7sin(2r) +c 


2 4 
was obtained by hand. The equivalence of this 


and the Mathcad expression follows from the 
double-angle formula sin(2x”) = 2sin x cos z. 


= “3 2 
(d) [ze = 5 ne 
See Activity 2.5(b) (iii). 


(e) [ue du = Zuexp(3u) — 5 exp(3u) +¢ 


+1)+e 


(f) [o In(5z) dx = $2° In(5a) — kx? +¢ 


1 
—— dt = aresin(44) +c 
e) / (-F st) 
(The function arcsin is represented in Mathcad 
by ‘asin’.) 
Solution 5.5 


The answers are presented in a form as close as 
possible to the Mathcad output (which by default is 
given to 3 decimal places). Each answer agrees with 
that obtained by hand in the main text. 


2 
(a) / e' dt = exp(2) — 1 = 6.389 
0 
See Activity 4.4(b). 
m/4 
(b) / (cos(5x) + 2sin(52)) dx = iQ? 4 2 


= 0.541 
See Exercise 4.1(a). 


o fs —= du=3 


See Exercise 4.1(b). (This answer is the 
Mathcad outcome from symbolic evaluation. ) 


(d) / e’sint dt = 4 exp(m) + 5 = 12.070 
0 * 


See Exercise 4.1(c). 
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The answers are given by Mathcad to 3 decimal 
places. 


A 
(a) / dt = 0880 
0 


1 
b) 4 V1 — 2? dx = 1.571 
il 
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Solution 4.1 


(a) The solution curves are described and sketched 


in the following table. 


(Xo, yo) Description 


Rough U-shape; 
steeper to right of 
minimum than to left. 


Similar, but with 
minimum to the left 
and higher. 


Similar, but with 
minimum still further 
to the left and higher. 


Straight line 
eee Mew 


Downward curve with 
gradient decreasing. 


(b) The solution curves show three distinct types of 
behaviour: 


(i) Any solution curve through a point above 
the line y = —x — 1 remains above that line. It 
decreases to a minimum and then increases. 


(ii) The line y = —z — 1 is itself the solution 
curve through any point on that line. 


(iii) Any solution curve through a point below 


the line y = —x — 1 remains below that line and 


decreases. 


The solution curves through the points (—3, —1), 


(—1,0) and (4,2) are of types (iii), (ii) and (i), 
respectively. 


Solution 4.2 


(a) (ii) The slope of the direction field at a point 
(xo, Yo) appears to depend on the choice of x9 
alone and not on that of yo. Correspondingly, 
there will be solution curves of just one type, 
with any two such curves differing only by a 
vertical translation. (This is a consequence of 
the fact that the function f(x,y) in this case 
depends only on 2.) 


(iii) All choices for zo and yo give a solution 
curve which increases and decreases alternately, 
but has a rising trend. Each such curve may be 
obtained from that shown below by a vertical 
translation. 


Figure S3.1 


(b) (ii) The effect of the direction field on solution 
curves in this case is not so clear-cut. However, 
it appears that there may be three types of 
solution curve, as detailed below. 


(iii) There are solution curves that cross the 


g-axis and are increasing, as in the graph below. 


If x9 = 0, then such curves are obtained by a 
choice of yo such that —1.25 < yo < 1.25. 


Figure S3.2 


SOLUTIONS TO ACTIVITIES 


There are solution curves that lie above the 
x-axis, each having a minimum in the second 
quadrant, as in the graph below. If x9 = 0, then 
such curves are obtained by a choice of yo such 
that aq > 1.26. 


Figure $3.3 


There are solution curves that lie below the 
x-axis, each having a maximum in the fourth 
quadrant, as in the graph below. If xo = 0, then 
such curves are obtained by a choice of yo such 
that yo < —1.26. 


Figure S3.4 


Solution 4.4 
(a) The graph extends to x = 9, but the portion to 


the left of x = 5 remains unchanged. In 
particular, the step size h does not alter, though 
the number N of steps does (so that the graph 
reaches = 9 rather than x = 5). 


A table of the values obtained is below. 


h ! 0-02-01 01 G00! 0.000) 
Un 3.016 3.001 3.002 2.887 2.783 2.772 2.771 


The values of yy appear to be converging and, 
given the level of agreement between the last two 
estimates, it seems reasonable to deduce from 
them an estimate for y(9) that is accurate to one 
decimal place, that is, y(9) = 2.8. (In fact, it 
looks likely that y(9) = 2.77 to 2 decimal places.) 
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Solution 4.5 


(a) 


The table below gives the values obtained using 
Euler’s method, from file 121C3-02. 


h 1 0.5 0.1 0.01 0.001 0.0001 
yn - 2.000-2.004 2.021 2.026 2027 207 


The agreement of the last two estimates suggests 
that y(6) = 2.027 to three decimal places. 


For the given function y = « — 4+ 4e'~*, we 
have 

yl) = 1 +44 42" = 1, 
so that the initial condition of the problem in 
part (a) is satisfied. 
The derivative of the given function is 

dy 


eS 1—2x 
dx Ertan 


whereas we have 
e—-y—-3=a2—-—(4%—4+4e'*)-3 
= 1—4e"*; 
Hence the given function also satisfies the 
differential equation from part (a), and therefore 


satisfies all the conditions of the initial-value 
problem. 


The value of this solution function at x = 6 is 
y(6) =6—4+4 4e' © = 24 de ~ 2.027. 


This value agrees to 3 decimal places with that 
obtained using Euler’s method. 
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